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ABSTRACT 

We compute two-point correlation functions and measure the shear signal due to 
galaxy-galaxy lensing for 80 000 optically identified and 5 700 radio-loud AGN from 
Data Release 4 (DR4) of the Sloan Digital Sky Survey. Halo occupation models are 
used to estimate halo masses and satellite fractions for these two types of AGN. The 
large sample size allows us to separate AGN according to the stellar mass of their host 
galaxies. We study how the halo masses of optical and radio AGN differ from those 
of the parent population at fixed M*. Halo masses deduced from clustering and from 
lensing agree satisfactorily. Radio AGN are found in more massive halos than optical 
AGN: in our samples their mean halo masses are 1.6 x 10 13 and 8 x lO 11 ^ 1 Mq, 
respectively. Optical AGN follow the same relation between stellar mass and halo 
mass as galaxies selected without regard to nuclear properties, but radio-loud AGN 
deviate significantly from this relation. The dark matter halos of radio-loud AGN are 
about twice as massive as those of control galaxies of the same stellar mass. This boost 
is independent of radio luminosity, and persists even when our analysis is restricted 
to field galaxies. The large-scale gaseous environment of the galaxy clearly plays a 
crucial role in producing observable radio emission. The dark matter halo masses that 
we derive for the AGN in our two samples are in good agreement with recent models 
in which feedback from radio AGN becomes dominant in halos where gas cools quasi- 
statically. 

Key words: galaxies: active - galaxies: haloes - galaxies: formation - gravitational 
lensing - dark matter -large-scale structure of Universe 



1 INTRODUCTION 

It is now widely accepted that galaxies form by the cool- 
ing and condensation of baryons within a m erging hierarchy 
of dark matter halos (|White fe Rees|[l97Sl ). Processes other 
than cooling also influence the relationship between galaxies 
and their halos. "Feedback," both from supernova explosions 
and from energy liberated during the accretion of material 
onto a central supermassive black hole, is currently under 
considerable scrutiny and debate, because it is believed to 
play a very important role in regulating the fraction of avail- 
able baryons that end up in galaxies. 

Over the past few years, considerable effort has been 
devoted to obtaining quantitative constraints on the rela- 
tionship between galaxies and their dark matter halos us- 
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ing a variety of different methods. Models that describe 
the evolved, non-linear dark matter distribution in terms 
of its halo building b l ocks (so-call e d halo models, e.g. 
Peacock fc Smith! |2000| : ISeliakl |2000|; iBerlind fc Weinberg! 
20021 : ICoorav fc Shethll2002l; lYang et al.ll2003l'l or direct N- 
simulations (e.g. 



body simulations (e.g. Kauffma nn et all] ] 
199Sl ; lKauffmann et aljl999l : iBenson et alJ li 



19971: Uing et al 
2000 ; I Yang et al 



2003), can be used in conjunction with the measured clus- 



tering amplitude of galaxies to constrain the relationship 
between the galaxies and their host halos. These constraints 
should be regarded as indirect, in part because of the need 
for a cosmology-dependent conversion from galaxy bias to 
halo mass. 

Weak lensing around galaxies (or galaxy-galaxy lens- 
ing, hereafter g-g lensing) provides a direct probe of the 
dark matter that surrounds galaxies (for a review, see 
iBartelmann fc S chneideJ l200ll ). Gravitational lensing in- 
duces tangential shear distortions of background galaxies 
around foreground galaxies, allowing direct measurement of 
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the galaxy-mass correlation function around galaxies. The 
individual distortions are small (of order 0.1 per cent), but 
by averaging over all foreground galaxies within a given sub- 
sample, we obtain high signal to noise in the shear as a 
function of angular separation from the galaxy. If we know 
the lens redshifts, the shear signal can be related to the pro- 
jected mass density as a function of proper distance from the 
galaxy. Thus we can observe the averaged DM distribution 
around any given galaxy sample. 

These techniques have been applied to study how 
galaxies with different properties such as luminosity, stel- 
lar mass, colour, spectral type and morphology populate 
dark matter halos of different masses dHoekstra et al.ll2005l ; 
iHevmans etaHbOOd iMandelbaum et al.ll2006h . The derived 
relations serve as important constraints on models of galaxy 
formation, but they do not directly constrain the physical 
processes that are responsible for creating the relations in 

the first place. 

For example, lYang et al.l (|2005l ) use the occupation 
statistics of a group catalogue to show that the mean lu- 
minosity of halo central galaxies scales with halo mass as 
L c oc M 2 / 3 for halos less massive than 10 13 K~ 1 Mq , and 
flattens to a much shallower relation (L c oc M 1 ^ 4 ) for more 
massive halos. It has been proposed that this characteristic 
scale reflects the imprint of feedback from radio- loud AGN, 
which heat the gas in massive halos and prevent it from cool- 
in g and conden s ing on to t he central galaxy. In the models 
of ICroton etafl (|2006l 1 and iBower et al.l l|20061 ), this "radio 
mode" feedback operates in halos with masses greater than 
~ 3 x 10 11 h~ x Mq, where cooling times are long compared 
to the free-fall time and gas co ols quasi-hydro s taticall y. On 
the other hand, the work of ISpringel et al.l l|2005al) and 
iHopkins et ail (|2005bl lah has focused on the role of optically 
luminous AGN in expelling gas from galaxies and regulat- 
ing the rate at which they are able to form stars. In these 
models, the major growth phases of black holes and the trig- 
gering of optically luminous AGN occur when two galaxies 
that contain sufficient cold gas merge with each other. The 
triggering thus does not depend directly on the mass of sur- 
rounding dark matter halo. 

In order to constrain the importance of these pro- 
cesses, it is important to understand how the AGN them- 
selves are related to the surrounding dark matter distri- 
bution. The large-scale clustering amplitude of luminous 
quasars has been accurately measured using tens of thou- 
sands of such ob jects drawn from th e 2dF and Sloan Dig- 
ital Sky Survey. ICroom et al.l (120051 ) use a sample based 
on the 2dF to conclude that quasars from 0.5 < z < 2.5 
inhabit dark matter halos with a characteristic mass of 
~ 3 x 10 12 hT 1 Mq, and that this m ass does not depend 
on redshift or on quasar luminosity. iMvers et all (|2007al ) 
use a photometrically-identified quasar sample from SDSS 
over a similar redshift range to estimate a typical mass of 
~ 5 x 10 12 /i -1 Mq, without a robust detection of luminosity- 
dependent bias at fixed redshift (for which there is only a 
marg inal detection in th e 2dF sample, Porciani fc Norberg 

20061 ; IMvers et al 



Seve ral studies dHennawi et al 



2007b, 2008) have also probed the very small-scale cluster- 



ing of quasars, using pairs of binary quasars, which are a 
probe of how the local environment affects quasar activity, 
though no clear co nsensus arises from these studies. The 
ICroom et al.l (|2005l ) results about halo masses were recently 



ex tended to redshifts z > 3 bylShen et al.l l|2007l ) and z < 0.6 
by iPadmanabhan et al.l (120081 ). The lack of evolution of the 
halo masses from the ~ 1O 12 /i -1 M0 scale is all the more 
remarkable considering that the nonlinear mass evolves by 
well over a factor of ten over the full redshift range probed 
by these studies. 

Clustering measurements of radio-loud AGN have been 
considerably less accurate because many fewer redshifts have 
been available. Studies of the angular clustering of radio 
sources drawn from the wide-area surveys such as the NRAO 
VLA Sky Survey (NVSS) or the Faint Images of the Radio 
Sky at Twenty-centimetres (FIRST) survey show that radio- 
loud AGN are considerably more strongly clustered than 
quasars. The est imated halo masses a re typically around 
10 13 - 10 14 M© (|Overzier et al.l l2003h . iMagliocchetti et all 
(2004) computed the redshift-space correlation function for 
820 nearby radio sources with redshifts from the 2dF, and 
derived a characteristic halo mass of 1O 13 4 M . No difference 
was found in the clustering properties of AGN with different 
radio luminosities. 

This paper focuses on a sample of 80 000 optically iden- 
tified AGN and 5 700 radio-loud AGN drawn from the Data 
Release 4 (DR4) of the Sloan Digital Sky Survey. We com- 
pute two-point correlation functions and measure the shear 
signal due to g-g lensing for these two samples. The large 
sample size allows us to split the AGN into different bins in 
stellar mass and study how the dark matter halo masses 
of optical and radio AGN differ at a fixed value of M, . 
Nearby radio-loud AGN have been shown to have signifi- 
cantly higher stellar masses than optic ally-identified AGN 
at the same redshift l|Best et al . 2005b), so it is important 
to understand whether the derived halo masses simply track 
this difference or whether the dark matter halo affects either 
the ability of an accreting black hole to produce a jet or the 
detectability of the jet at radio wavelengths. 

We also create control samples of non-AGN that are 
matched in stellar mass, redshift and morphology, and we 
use these control samples to investigate whether the halo 
masses of active galaxies differ from those of their counter- 
parts chosen irrespective of their level of nuclear activity. 
Finally, the fact that the clustering and weak lensing anal- 
yses are carried out on the same set of galaxies allows us to 
evaluate the consistency of our constraints on dark matter 
halo mass obtained using the two methods. 

We begin by outlining the theory behind the lensing 
and clustering measurements in section [21 We then describe 
the data used for the analysis, and the analysis procedure, 
in section [3] The results for optical and radio-loud AGN 
are presented in section 3] including lensing and clustering 
separately, followed by a joint analysis. We then summarize 
the key results and discuss their implications in sections [S] 
and [51 respectively. 



2 THEORY 

2.1 Galaxy-galaxy lensing 

Galaxy-galaxy weak lensing provides a simple way to probe 
the connection between galaxies and matter via their cross- 
correlation function 



£,gm{r) = {5 g (x)5 m (x + r)) 



(1) 
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where 5 g and S m are overdensities of galaxies and matter, 
respectively, and in practice the mean is taken over some 
survey volume (in theory it is the average over a whole dis- 
tribution, but we can only estimate its value using the fixed 
volumes that are available in reality). We will interchangably 
express correlation functions as functions of vectors (r) and 
scalars (r) because of the assumption of statistical isotropy. 

This cross-correlation can be related to the projected 
surface density 



(0P 



(for r 2 — R 2 +\ 2 ), where we ignore the radial window, which 
is much broader than the typical extent of the lens. This 
surface density is then related to the observable quantity 
for lensing, 



AE(-R) = 7t(-R)£ c = £(< R) - E(7?), 



(3) 



where the second relation is true only for a matter distri- 
bution that is axisymmetric along the line of sight. This 
observable quantity can be expressed as the product of two 
factors, a tangential shear 74 and a geometric factor, the 
critical surface density 

E^A ^^f. „ (4) 



4ttG D l D ls (1 + z l ) 2 

where Dl and Ds are angular diameter distances to the lens 
and source, Dls is the angular diameter distance between 
the lens and source, and the factor of (1 + zl)~ 2 arises due 
to our use of comoving coordinates. For a given lens redshift, 
E^ 1 rises from zero at z s — zl to an asymptotic value at 
z s S> zl; that asymptotic value is an increasing function of 
lens redshift. 

In practice, we measure the g-g weak lensing signal 
around a stacked sample of lenses to obtain the average 
AE(.R) for the whole sample. This stacked lensing signal 
can be split into two terms that dominate on different scales. 
The 1-halo or Poisson term, which is determined by the dark 
matter halo in which the galaxy lives, dominates on scales 
typically below ~ lh~ Mpc. The halo-halo term, which is 
determined by correlations between the galaxy and other 
dark matter halos, dominates on larger scales. The 1-halo 
term can be further split into two contributions. For central 
galaxies, which reside at the peak density of a dark matter 
halo that is not contained within another halo (a host halo) , 
the Poisson term is simply determined by the matter density 
of that host halo, p(r). For satellite galaxies, which reside 
in dark matter subhalos, there is a contribution from the 
density of the dark matter subhalo, but there is also a term 
on hundreds of kiloparsec scales due to the cross-correlation 
between the galaxy position and the host dark matter halo. 
Consequently, the lensing signal on <~ 0.3/i -1 Mpc scales 
tells us about the dark matter halo in which the galaxy re- 
sides; the signal from ~ 0.3 - lh~ Mpc reveals the local 
environment of the galaxy; and the signal on larger scales 
indicates the large-scale correlations of the galaxy sample. 

We interpret the lensing signal statistically using a halo 
model, which allows us to determine both the typical halo 
mass M cen t for central galaxies in our galaxy sample, and 
also the satellite fraction a (the fraction of the sample lo- 
cated in subhalos within some more massive host dark mat- 
ter halo). In this simple formulation of the halo model, we 



assume that all central galaxies in our sample have a sin- 
gle halo mass M cen t, and that all satellite galaxies are dis- 
tributed in halos with M > 3M cen t with the number in a 
halo of given mass above this threshold oc M. Tests of this 
halo model f ormulation using the lensin g signal from N-body 
simulations (|Mandelbaum et a l. 2005b) clearly indicate that 
the best-fitting M cen t and a recover the true values to within 
~ 10 per cent, provided that the distribution of central halo 
masses is relatively narrow (FWHM typically a factor of 6 or 
less). For a broader distribution of central halo masses, the 
best-fitting M cen t lies between the median and the mean of 
the distribution, but may differ from either one by as much 
as a factor of two. For more details of t his hal o model and its 
assumptions, see iMandelbaum et all <|2005bh . In this work, 
we use this halo model without applying any correction for 
the (unknown) scatter between galaxy stellar mass and dark 
matter halo mass, but we will discuss the extent to which 
our assumptions about the width of the central halo mass 
distribution are likely to be correct. 



2.2 Galaxy clustering 

The clustering of galaxies is usually q uantified using the two- 
point correlation function (2PCF, e.g.|Pecbles 1980), defined 
by 



dPn = n 2 [l + £(r)]dVidV 2 . 



(5) 



Here n is the mean number density of galaxies, and dVi 
and dV2 are the volumes of two infinitesimally small spheres 
centered at x\ and X2 with distance of r = — 3\. By defi- 
nition, dP\2 is the joint probability that a galaxy lies in each 
of the spheres, and so the 2PCF £(r) represents the excess 
probability of finding two galaxies separated by a distance r , 
compared with the result obtained for a uniform random dis- 
tribution. If £(r) > 0, then galaxies are said to be clustered. 
In galaxy redshift surveys, the 2PCF is measured in redshift 
space and usually expressed as functions of separations per- 
pendicular (r p ) and parallel (n) to the line of sight. In many 
cases, the projected two-point correlation function, w p (r p ), 
is the more useful quantity, because it does not suffer from 
redshift-space distortions, and is thus directly related to the 
real-space correlation function. The 2PCF is also simple to 
compute and can be easily compared with the predictions 
of theoretical models. 

The amplitude of the correlation function on scales 
larger than a few Mpc provides a direct measure of the mass 
of the dark matter halos that host the g alaxies t h rough the 
halo mass - bias relation. As shown in iLi et all (|2008al lbh . 
the amplitude of the correlation function on scales < 100 
kpc can serve as a probe of physical processes such as merg- 
ers and interactions. On intermediate scales, the shape of 
the correlation function is sensitive to how galaxies are dis- 
tributed within their dark matter halos. 

The clustering signal must also be interpreted using 
some form of halo model. Following the approach adopted in 
our previ ous work , we int erpr et clustering result s using the 
models of lLi et all l|2006al ) and I Wang et all (|2006l ) which are 
based on direct N-body simulations. We have constructed a 
set of 10 mock galaxy catalog ues from the Millennium Sim- 
ulation (jSpringel et alll2005bT ) with exactly the same obser- 
vational selection effects as the SDSS DR4. The Millennium 
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Simulation uses 10 10 particles to follow the dark matter dis- 
tribution in a cubic region 500/i _1 Mpc on a side. The cos- 
mological parameters assumed are fi m = 0.25, S7a = 0.75, 
erg = 0.9 and h = 0.73. We adopted the positi ons and veloc- 
ities o f the galaxies given in the catalogue of ICroton et al.l 
(2006), who implemented a semi-analytic model in order to 
track the formation and evolution of galaxies in the sim- 
ulation. Physical properties of the galaxies, such as stel- 
lar masses and AGN status, are not taken from the semi- 
analytic model, however, but instead are assigned to each 
model galaxy using parametrized functions. The main such 
function relates the stellar mass of the galaxy to the mass of 
the halo at the epoch when the galaxy was last the central 
dominant object in its own halo, including scatter in that 
relation. Tests have shown that this procedure allows us to 
match accurately both the stellar mass function of SDSS 
galaxies and the shape and amplitude of their two-point 
correlations as a function of stellar mass (ILi et al . 2006a; 
IWang et alj|2006h . 



In I Li et al] (|2006bl . hereafter L06), we adapted this 
halo model to interpret the clustering of optically-identified 
AGN. In that paper, we computed the correlation functions 
of AGN and of control samples of inactive galaxies that had 
the same redshift and stellar mass distribution as the ac- 
tive galaxies. We found that on scales between 100 kpc and 
1 Mpc, AGN are clustered more weakly than the inactive 
sample. We then introduced a simple model in which the 
probability of a galaxy of given stellar mass to be an AGN 
is enhanced if it is the central galaxy of its own halo, and 
showed that this model could provide a good fit to the data. 
In the best-fitting model, 84 per cent of all optical AGN are 
located at the centres of their own dark matter halos (i.e., 
fcen = 0.84), whereas this is true for only 73 per cent of 
inactive galaxies. 



We emphasize that while in principle, it would be easy 
to assume that a derived from the lensing analysis is sim- 
ply 1 — fcen from the clustering analysis, the relationship is 
not as straightforward as this. The clustering analysis starts 
from the assumption that the AGN and control samples de- 
rive from a parent population which matches the statisti- 
cal properties of galaxies as a function of stellar mass (i.e., 
the stellar mass function and the clustering as a function of 
mass) . Any halo model parameters such as f cen are then in- 
troduced as a way of matching the observed AGN clustering 
by modifying the probability that each galaxy in the parent 
population is an AGN. In contrast, the lensing analysis does 
not assume that the AGN and control samples stem from 
identical parent populations. In future analyses with larger 
datasets, we will jointly model clustering and lensing with 
the same assumptions at the outset; here, we simply use 
pre-existing analysis pipelines that output quantities that 
should be reasonably (but not exactly) comparable. Even in 
the absence of a completely unified approach to modeling, 
there are many valuable conclusions that can be drawn from 
the lensing and clustering signals (e.g., if both the clustering 
and lensing signals for a particular sample are comparable to 
the signals for the controls, or if they are both quite different 
from the signals for the controls). 



3 DATA AND SIGNAL MEASURES 
3.1 Overview of SDSS 

The d ata used here are obtained from the SDSS (|York et al.l 
2000), an ongoing survey to image roughly n steradi- 
ans of the sky, and follow up approxima tely one million 
of the detected objects spectroscopically (lEisenstein et al.l 
l200ll ; iRichards et al.l I2OO2I : IStrauss et all I2OO2I ) . The imag- 
ing is carried out by drift-scanning the sky in ph oto- 



metric conditions dHogg et al.l 200 ll; Ivezic et al. 2004), in 



five bands (ugriz) (|Fukugita et al.lll996l ; ISmith et al.ll2002l) 



using a specially-designed wide-field camera ( Gunn et alj 
1998). These imaging data are used to create the source 
catalogue that we use in this paper. In additi on, objects are 
target ed for spectroscopy using these data (|Blanton et al.l 
2003oj) and are observe d with a 640-fibe r spectrograph 
on the same telescope l|Gunn et al.1 [2006). All of these 
data are processed by completely automated pipelines 
that detect and measure photometric properties of ob- 



200 ll iPier et al.1 120031; iTucker et al. 


2006). The SDSS has 


had seven major data releases ( 


Stoughton ct al. 


2002; 


Abazaiian et al. 2003L 2004 2005; 


Finkbcincr ct al. 


2004; 


Adclman-McCarthv et all 20061. 2007, 20081). In this Da- 



per we use data from the fou rth of these releases(DR4; 
lAdelman-McCarthv et~ai1l2006l ). 



3.2 The AGN and control samples 

3.2.1 Optically-identified AGN 

The sample of opti cally-identified A GN is the same as that 
analyzed in L06 and lLi et all l|2008bl ). in which the clustering 
of AGN on a variety of different scales was studied. 

The base sample is composed of ~ 4 x 10 5 objects for 
which data are publicly available through DR4 and which 
have been spectroscopically confirmed as galaxies with r- 
band magnitudes in the range 14.5 < r < 17.6, redshifts 
in the range 0.01 < z < 0.3, and absolute magnitudes in 
the range —23 < Mo.i r < —17. Here r is the r-band Pet- 
rosian apparent magnitude corrected for foreground extinc- 
tion, and Mo.i r is the r-band absolute magni tude corrected 
to its z 



0.1 value using the fc-correct code of lBlanton et al 
(2003a) and the luminosity evolution model of lBlanton et al 
(|2003bl ). A sample of ~ 80,000 AGN are selected from the 
subset of these galaxies with S/N > 3 in the four emission 
lines [O m]A5007, H/3, [N ii]A6583 and Ha , following the 
criteria proposed bv lKauffmann et al.l (2003). In the follow- 
ing analysis, we occasionally divide our sample into "weak" 
and "powerful" AGN using the quantity L[OIII]/Mt,h, where 
L[OIII] is the extinction-corrected [OIII] line luminosity of 
the AGN and Mth is the black hole mass estimated from 
th e velocity dispers i on of the galaxy using t he relation given 
in iTremaine et ail ((2002) . As discussed in iHeckman et al.l 
l|2004l ). this quantity can be viewed as a measure of the ac- 
cretion rate onto the black hole relative to the Eddington 
rate. 

We have also constructed two sets of 20 different control 
samples from the full parent sample of galaxies by matching 
a number of physical parameters regardless of nuclear activ- 
ity. For the first set, four physical parameters are matched: 
redshift (z), stellar mass (M„), concentration (R90/R50), 
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and stellar velocity dispersion (it*). For the second set, the 
4000 A break strength (D4000) is also matched. The match- 
ing tolerances are Acz < 500 km s~\ AlogM, < 0.L 
Aa* < 20 km s _1 , AR 90 /R 50 < 0.1 and AD4000 < 0.050 
In the first case, 28 per cent of the control galaxies are also 
included in the AGN sample; in the second case, the over- 
lap fraction is higher, 37 per cent. In both cases, the match 
fraction is a slightly increasing function of stellar mass. 



3.2.2 Radio-loud AGN 

The National Radio Astronomy Obser vatory (NRAO) Very 
Large Array (VLA) Sky Survey fNVSS; ICondon et al.|[l998r ) 
and the Faint Images of the Radio Sky a t Twenty centime- 
ters (FIRST) survey l|Becker et al.l fi9951 are two radio sur- 
veys that have been carried out in recent years using the 
VLA radio synthe sis telescope at a frequency of 1.4 GHz. 
iBest et al l (l2005bh identified radio-emitting galaxies within 
the main spectroscopic sample of the SDSS data release 2 
(DR2) by comparing these galaxies with a combination of 
these two surveys. The use of two radio surveys allowed a 
radio sample to be constructed that was both reasonably 
complete (~95 per cent) and highly reliable (it is estimated 
that ~99 per cent of the sources in the catalogue are genuine 
radio galaxies rather than false matches). In this paper, we 
use an updated catalogue of 5 712 radio galaxies based on 
the SDSS data release 4 (DR4). The sample spans the same 
redshift range as the sample of optically-selected AGN, and 
the radio luminosities of the AGN range from 10 23 to 10 26 
W Hz -1 (i.e., they are mainly FRI type systems). 

We have used the parent galaxy catalogue to construct 
a set of five control samples that are closely matched in 
redshift, stellar mass and stellar velocity dispersion. The 
matching tolerances are the same as used to construct the 
optical AGN control samples. We do not additionally match 
in D4000, because almost all galaxies at the relevant stellar 
masses are red, with a strong 4000 A break. Of the control 
sample, 10.3 per cent of the galaxies are radio-loud AGN; for 
our higher stellar mass subsample (log(M*/MQ) ^ 11.44), 
the radio- loud AGN fraction is 16 per cent, whereas for the 
lower M, subsample, it is 7 per cent. 



3.3 Lensing analysis 

3.3.1 Lensing source ca 

The source sample used for th e lensing analysis is t he sam e 
as that originally described in iMandelbaum et al.l (l2005aT ) . 
This source sample includes over 30 million galaxies from the 
SDSS imaging data with r-band model magnitude brighter 
than 21.8, with shape measurements obtained using the 
REGLENS pipeli ne, including PSF cor rection done via re- 
Gaussianization (|Hirata fc Seliakl [2003) and with cuts de- 
signed to avoid various shear calibration biases. In addition 
to these, there are also uncertainties due to photometric 
redshifts and / or redshift distributions of background galax- 
ies, as well as due to other issues affecting the calibration 



1 This procedure is identical to that in L06, except that the con- 
trol galaxies are selected from the full parent sample, rather than 
from the subset of inactive galaxies. 



of the lensing signal, such as the sky subtraction uncer- 
tainties, intrinsic alignments, magnification bias, star-galaxy 
separation, and seeing-dependent systematics. The overall 
calibration uncertainty wa s estimated to be eight per cent 
l|Mandelbaum et al.ll2005ah . though the redshift calibration 
component of this systematic error budget has recently been 
decreased due to the avail ability of more spectroscopic data 
ijMandelbaum et al.ll2008T ). With a total estimated lensing 
calibration uncertainty of ~ 5 per cent, this systematic is 
subdominant compared to the statistical error and to the 
uncertainty derived from the model used to interpret the 
lensing signal. 



3.3.2 Lensing signal computation 

Here we briefly describ e the computation of the le nsing sig- 
nal; for more detail, see lMandelbaum et al.l l|2005al ). For each 
lens, we identify sources within 46 logarithmically-spaced 
annuli around the lens (in comoving transverse separation) 
from 20 fe _1 kpc to 2 /i _1 Mpc. The tangential ellipticity of 
the source relative to the lens is measured, in order to es- 
timate the tangential shear. Lens-source pairs are assigned 
weights according to the error on the shape measurement 



Wis 



al + a 2 SN 



(6) 



where a% N is the intrinsic shape noise and a s is the mea- 
surement error on the source galaxy ellipticity. The factor of 
E~ 2 optimally weights the signal by the noise in AE rather 
than in the shear. 

Once we have computed these weights, we compute the 
lensing signal in each radial bin as a summation over lens- 
source pairs via: 



AE(R) 



Y,is w 't 



(7) 



where the factor of 2 and the shear responsivity 1Z relate our 
definition of ellipticity to the shear, using the formalism in 
iBernstein fc Jarvisl (|2002t ). In practice, 1Z ~ 1 — e 2 ms w 0.86. 

There are several additional procedures that must 
be done when computing th e signal (for more detail, 
see IMandelbaum et al.ll2005al h First, the signal computed 
around random points must be subtracted from the signal 
around real lenses to eliminate contributions from system- 
atic shear. In practice, this correction is negligible for the 
scales used in this work. Second, the signal must be boosted, 
i.e. multiplied by B(R) — n(R)/n ra nd(R), the ratio of the 
number density of sources relative to the number around 
random points, in order to account for dilution by sources 
that are physically associated with lenses, and therefore not 
lensed. 

To determine errors on the lensing signal, we divide 
the survey area into 200 bootstrap subregions, and gen- 
erate 2500 bootstrap-resampled datasets. These bootstrap- 
resampled datasets are also crucial for determining the sta- 
tistical significance of differences between correlated sub- 
samples of galaxies, because fitting the signal to the halo 
model on each resampled dataset allows us to determine 
how much any overlap between two galaxy samples leads to 
a correlation between the best-fitting halo model parameters 
for the two samples. 
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Figure 1. The galaxy-galaxy lensing signal for optical AGN (top) and control galaxies (bottom) as a function of transverse separation. 
Points are the data, and lines are the best-fitting halo models. The left panels show the sample split by stellar mass as follows: logM* < 
10.6 (black solid), 10.6 ^ logA/» < 11 (red dotted), logAf, ^ 11 (blue dashed). The right panels show the sample split by L[Olll}/M bh 
into the lower half of the sample (black solid) and upper half (red dotted). 



The lensing signal is presented in comoving coordinates, 
with angular diameter distances computed assuming a flat 
ACDM universe with Q m = 0.3 and f2 A = 0.7. The halo 
model used to interpret the lensing signal assumes erg = 0.9. 
In the units used, Ho scales out of everything, so our results 
are independent of this quantity. The central halo mass def- 
inition for this paper is the mass within which the spherical 
overdensity is 200p crit , which is roughly 35 per cent lower 
than the mass definition used fo r previous lensing analyses 
using this halo model formalism (|Mandelbaum et al.ll2005bl . 
2006). This change in halo mass definition was made to 
match the mass definition for the clustering analysis. 



3.4 Clustering analysis 

3.4-1 The reference galaxy sample 

In this paper, the clustering of AGN (or control galaxies) 
is quantified by the projected two-point cross-correlation 
function (2PCCF), w p (r p ), which is estimated by cross- 
correlating the AGN (or control) samples described above 
with a reference sample of galaxiesff] The reference galax- 
ies are selected from sample dr4 of the New York Uni- 
versity Value Added Galaxy Catalogue (NYU-VAGC), 
which is based on SDSS DR4, publicly available at 

2 We use the notation r p for the transverse separation in the 
clustering analysis, and the notation R for the same quantity in 
the lensing analysis. The main reason is to maintain notational 
consistency within previous work. 
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http://sdss.physics.nyu.edu/vagc/, and is described in de- 
tail in lBlanton et al.l l|20051 ). The reference sample contains 
292,782 objects that are identified as galaxies from the Main 
sample and have 0.01 ^ z ^ 0.3, 14.5 < r < 17.6 and 
—23 < Mo.i r < —17. This sample has formed the basis of 
our recent investigations of the clustering properties of dif- 
ferent classes of galaxies (L06, iLi et al.|[2007al . l2008al lbh . 



3.4-2 Clustering measures 

Our methodology for computing correlation functions has 
also been described in detail in our previous papers. Random 
samples are constructed with the same selection function 
as the reference sample, as described in detail in ILi et al.l 
(2006a) (but note the slight differences mentioned here in 
%T2T) . The redshift-space 2PCCF £(r p ,ir) between AGN 
(or control galaxies) and the reference sample is then calcu- 
lated using the estimator presented in L06, 



Trends with stellar mass for optical AGN 



N R QD{r v ,n) 
^ Tp ' n) ~ N D QR(r p ,n) L ' 



(8) 



where r p and tt are the separations perpendicular and par- 
allel to the line of sight; Nd and Nr are the number of 
galaxies in the reference sample and in the random sample, 
with Nr/Nd = 10 throughout this paper; QD(r p ,ir) and 
QR(r p ,n) are the cross pair counts between AGN/control 
and the reference sample, and between AGN/control and 
the random sample, respectively. Finally, the redshift-space 
projected 2PCCF w p (r p ) is estimated by integrating £(r p , it) 
along the line-of-sight direction: 



w p (r p ) 



t(r p ,Tr)dn 



^(r p ,TTi)ATVi- 



(9) 



Here Tr max — 40/i _1 Mpc, and the summation for comput- 
ing w p (r p ) runs from 7Ti = —39.5 h _1 Mpc to 7rso = 39.5 
h -1 Mpc, with AiVi = 1 h _1 Mpc. We have also corrected 
carefully for the effect of fibre collisions; a description and 
tests of the method can be found in L06. As will be de- 
scribed in more detail in § 14.2.11 error estimates come from 
the variance in w p (r p ) between 100 mock catalogues. 

The clustering computation assumes the same flat 
ACDM universe with fl m = 0.3, Qa = 0.7 and erg = 0.9 as 
for the lensing analysis. Our results are presented in units 
of ir 1 Mpc with h = l. 



4 RESULTS 

4.1 Optical AGN 

Results and interpretation for the galaxy clustering signal of 
the optical AGN have been presented in L06, and are also 
briefly described in £12.21 Consequently, here we present only 
the galaxy-galaxy lensing signal and its interpretation for 
this sample. In section [5j we compare the halo masses and 
satellite fractions estimated through lensing with the same 
quantities estimated through clustering. 

In Fig. [1] we show the g-g lensing signal for the optical 
AGN sample split by stellar mass and by the accretion rate 
per unit black hole mass (L[OIII] / Mbh) ■ Results are shown 
for both the AGN and the control samples; results for the 
control samples with the same distribution of D4000 are not 
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Figure 2. The best-fitting central halo masses M cen t (top) and 
satellite fractions a (bottom) as a function of stellar mass for the 
optical AGN and the two control samples as labelled on the plot. 



shown, because they are statistically consistent with the re- 
sults for the control samples where D4000 is not matched. 

There are several clear trends in this figure. First, the 
g-g lensing signal on small scales (< 0.3h~ Mpc) shows that 
the halo mass for central galaxies increases with stellar mass 
and decreases with L[OIII]/Mi,fe. This conclusion is true for 
both the AGN and for the control samples. Second, the g-g 
lensing signal for AGN and controls in a particular sub- 
sample is quite similar; any differences are not statistically 
significant. 

We now consider the halo model interpretation of these 
results, represented by the best-fitting central halo mass and 
satellite fraction for each sample. These quantities are plot- 
ted for the optical AGN and the two control samples in Fig. [2] 
as a function of stellar mass. For reference, they are tabu- 
lated for all optical AGN and control subsamples, including 
the splits by L[Olll]/M bh , in Table Q] 

A few trends are evident from the plot and table. First, 
for the samples split by stellar mass, the differences in halo 
model parameters for the AGN and control samples are 
not statistically significant, as expected from Fig. [T] How- 
ever, the central halo mass shows a strong trend with stel- 
lar mass, consistent with the l ensing results for the gene ral 
galaxy population discussed in lMandelbaum et al.l l|2006l ). If 
we compare against the results in that paper after account- 
ing for the different halo mass definitions, we conclude that 
for the lower and middle stellar mass bins, the best-fitting 
central halo mass is consistent (within the noise) with the 
results for both early and late type galaxies, which have 
similar mean halo masses below stellar masses ~ 10 n MQ. 
For the highest stellar mass bin, our best-fitting central halo 
mass is more consistent with the results for late- type galaxies 
(lower by a factor of a few than that for early-type galax- 
ies). This result is consistent with the general tendency of 
these narrow-line AGN to be associated with galaxies with 
ongoing star formation. 
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Table 1. Best- fitting halo model parameters for fits to the optical AGN g-g weak lensing signal, with 68 per cent CL errors. 

Optical AGN Controls Controls with D4000 

Sample M cen t ol M cen t a M c 



IO^/i^Mq lO 12 ft -1 Af 10 12 h- 1 M, 







logM*/M Q < 10.6, (M*/(10 10 M Q )) = 2.3 0.13±0.08 0.26±0.04 0.31 ±0.16 0.27±0.05 0.29±0.15 0.33 ± 0.05 

10.6 < log M*/M© < 11, (M»/(10 10 M Q )> = 6.5 0.86 ± 0.21 0.18±0.04 0.77±0.20 0.22 ± 0.05 0.79 ± 0.21 0.19±0.05 

logM»/M Q > 11, (Af„/(10 10 M Q )) = 15.4 2.3 ± 0.4 0.19±0.04 3.0±0.5 0.20±0.04 2.8 ± 0.6 0.21 ± 0.05 

Lower half, L[OIII]/M 6h 1.1 ±0.2 0.25 ± 0.03 1.5 ± 0.2 0.25 ± 0.03 1.1 ±0.2 0.30 ± 0.03 

Upper half, L[OUl]/M bh 0.47 ±0.13 0.16 ±0.03 0.39 ±0.14 0.22 ± 0.03 0.52 ±0.16 0.17 ±0.03 




10 10 10" 10 12 10 13 10 14 10 15 



M halo [h"'M sun ] 

Figure 3. The distributions of central halo mass in the mock cat- 
alogues that are able to reproduce the clustering signal in lLi et all 
(2006b) in our three stellar mass bins. 

The satellite fractions decrease slightly from the low- 
est to middle stellar mass bin. Consistent with the results 
from L06, we find slightly lower satellite fractions for the 
optical AGN than for the control samples. Unlike for the 
galaxy clustering signal, this difference is not statistically 
significant. 

There is clearly a significant difference in the mean cen- 
tral halo mass and satellite fraction for the samples split 
at the median value of L[OIU]/Mbh- However, this quan- 
tity is itself correlated with stellar mass, so some of the 
trend derives from that correlation. For the lower half of 
the sample in OIII luminosity, the mean stellar mass is 
(M„> = 9.1 x 10 10 A'/ Q ; for the upper half, it is 7.3 x 10 10 M©. 
The results for the control samples with the same stellar 
mass distribution suggest that the difference in best-fitting 
central halo masses can be explained solely by this difference 
in stellar mass distributions. 

Finally, we discuss the broadness of the central halo 
mass distribution. As we have already noted, a broad central 
halo mass distribution would lead to the best-fitting central 
halo masses being an overestimate of the median mass, and 
underestimate of the mean mass. To assess whether this may 
be the case, we use the halo occupation models of L06 for 
optical AGN (see $2.2$ which can be used to derive a cen- 



tral halo mass distribution. This plot is shown in Fig. [3] 
As shown, the FWHM of the distribution for the two lower 
stellar mass bins is within the factor of ~ 6 needed for ac- 
curate estimation of the mean central halo mass. However, 
the distribution is sufficiently broad for the highest stellar 
mass bin that our estimate from the lensing signal is likely 
an underestimate of the mean, possibly by as much as 50 per 
cent. We do not apply a correction to determine the mean 
central halo mass for either the AGN or the controls in this 
bin, since there is significant systematic uncertainty in the 
correction factor itself. 



4.2 Radio AGN 

4-2.1 Galaxy Clustering 

We have measured the projected 2PCCF w p (r p ) of the radio 
AGN with respect to the reference galaxies, and compared 
this to the average result of the five control samples. The 
results are shown in Fig. [3] for the whole sample (circles in 
the left panel) and for two different ranges in stellar mass 
M». 

Radio AGN are more strongly clustered than control 
galaxies on all scales. The difference in w p (r p ) amplitude is 
a constant factor on scales smaller than ~ 1 Mpc, then rises 
slightly on larger scales. Fig.|3]also shows that the clustering 
amplitude of radio AGN increases with the stellar mass of 
the host galaxy. This result is consistent with the fact that 
more m assive galaxies are more strongly clustered (|Li et al.l 
l2006al ). We have also tested whether there is a dependence 
of the clustering amplitude on the radio luminosity of the 
AGN, and we find that at fixed M*, there is no significant 
effect. This finding is c onsistent with the recent results of 
iKauffmann et al.l {2008), who show that radio- loud AGN 
are in denser environments than control radio-quiet galax- 
ies, but that there is no dependence of local density on the 
radio luminosity of the AGN. 

We now use our mock catalogues (see Sj3]) to model the 
observed clustering measurements of radio AGN. We first 
tried to vary the fraction of radio AGN assigned to central 
versus satellite galaxies (as was done for the optical AGN). 
We found that if radio AGN are preferentially found in satel- 
lite galaxies, we can fit the data on scales smaller than a few 
Mpc, but the model then underpredicts the clustering am- 
plitude on larger scales. As we have discussed, the amplitude 
of the correlation function on scales larger than a few Mpc 
provides a direct measure of the mass of the dark matter 
halos hosting the ra dio AGN. Motiva t ed by the models of 
ICroton et all (|2006h and iBower et al.l (2006), we impose a 
lower threshold in halo mass, M™ m , as a second free pa- 
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Figure 4. Projected cross-correlation function w p (r p ) for radio AGN (green filled circles) in different stellar mass ranges compared 
to results for control samples selected without regard to AGN properties (black open symbols). Results for the best- fitting model are 
indicated as red (AGN) and blue (controls) shaded regions, where the width of the shaded regions corresponds to the 1 — a variance 
between 200 mock catalogues. Errorbars are significantly correlated (> fO per cent) between radial bins above ~ \h~ x Mpc. 



rameter of the model. In other words, radio-loud AGN are 
only found in dark matter halos more massive than M™ ln . 
The probability of a galaxy to be a radio-loud AGN depends 
not only on whether it is a central or satellite system, but 
also on the mass of its dark matter halo. However, the prob- 
ability that a particular galaxy is an AGN does not depend 
on the AGN status of its neighbors. While this step function 
in mass is undoubtedly an over-simplification, we adopt it 
as a first attempt at modeling to see if it is close enough to 
reality that the observations can be modeled in this way. 

We have generated a grid of 322 models by varying the 
two parameters, / cen and M™ n , with f cen ranging from 
0.40 to 0.84 with a step size of 0.02, and \og(M^ in /h' 1 M e ) 
ranging from 10.0 to 13.25 with a step size of 0.25. We have 
constructed 200 mock catalogues of radio AGN for each 
of the models. We measure w p (r p ) and its covariance ma- 
trix C at each grid point, and we compare the measure- 
ments to the SDSS results. In order to identify the best-fit 
model, we first calculate the likelihood of each parameter 
set, L(f cen , M^"), by 



L(f cen ,Mr n ) 

where 



(27T)" 



x/2 



det[C(/ cen ,M™ 



-0.5x^(/ce„,Af™ 1 ") 



(10) 



X 2 (/ ce n,Mr n )=X T C- 1 X. 

Here X = {Xj} (j = 1, m) is an m x 1 vector with 



1 \ - 



1JJp,i J ) model I Wp ) £ 



(11) 



(12) 



where n — 200 is the number of mock catalogues, m is 
the number of radial bins over which w p (r p ) is measured, 
w p{ r P,j)sDSS is the clustering amplitude at the j th radial 
bin as measured from the SDSS, and w Pt i(r Pt j) mo dei is the 
result at the j radial bin as measured with the i th mock 



catalogue. The m x m matrix C = {C'ij} = l,...,m) 
is the covariance matrix of the measurements from the 200 
mock catalogues, given by 



Ci.j — 



n - 1 



(13) 



where 



Yk,i — W Pt k(r p ,i)model (14) 

is the measurement at the i th radial bin from the k th mock 
catalogue, and 



(15) 



(Yi) = — > J Wp,k(r p ,i)model 
fc=l 

is the mean measurement at the i th radial bin over all mock 
catalogues. 

We define the best-fit model to be the one giving a 
minimum A computed as follows: 

A(f cen ,Mr n ) = -21nL(/ ceni Mr") (16) 
= X 2 (/cen, Mr™) + 2 ln{det[C(/ cen , Mr")]} + mln(2^). 

Note that this maximum likelihood estimate differs from the 
simple minimum of \ 2 if the determinant of the covariance 
matrix det[C] varies with the parameters. This variation is 
demonstrated in Figure [5] where we plot ln[det(C)] in the 
grid of the two model parameters, f cen and M™ 8 ™. As can 
be seen, the determinant of the covariance matrix does vary 
systematically from model to model. This variation is due to 
the fact that the covariance depends on the two- and four- 
point functions of the galaxy distributions, which clearly 
differ across the grid due to the different ways the halos are 
populated with AGN. However, the variation is relatively 
smooth, indicating that we have used enough mock cata- 
logues at each grid point (200) to determine the covariance 
matrix with a sufficiently small noise level. 
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Figure 5. Determinant of the covariance matrix of w p (r p ) (see Eg, I13|l . on the grid of the two model parameters, / cen and Af™ ln . f C en 
is the fraction of radio AGN that are hosted by the central galaxy of their own dark matter halo, and M™ n is the minimum mass for 
the halos that can host radio AGN. The contour levels are indicated at the right-hand side. 



We compare the measured w p (r p ) with the models in 
four radial bins centered at r p = 0.21, 0.65, 2.1 and 6.5 
/i -1 Mpc, with a step size of Alogr p = 0.5 (larger than the 
radial bins shown in the plots). The choice of 4 radial bins 
was motivated by tests showing that the covariance matri- 
ces evaluated from the mock catalogues were well-behaved 
in this case, whereas using a significantly larger number of 
radial bins causes the covariance matrices (a) to be noisier, 
and (b) to have peculiar patterns of correlations between 
bins suggestive of edge effects (in particular, strong corre- 
lations between certain adjacent pairs of radial bins that 
are not representative of the overall pattern of correlations). 
These particular radial bins were chosen to sample separate 
parts of the HOD, namely the central 1-halo term, the satel- 
lite 1-halo term, the transition between the 1- and 2-halo 
terms, and the 2-halo term (respectively). The minimum ra- 



dius was chosen because the data at smaller separations are 
rather noisy. 

Fig. [6] plots the contours of AA = A — A m ; n in the grid 
of the two parameters, when using the w p (r p ) measurements 
for the full radio- loud AGN sample (left panel of Fig.Q}. The 
1,2, and 3a confidence regions, computed for m — 4 and 2 
parameters, are indicated using solid, dashed and dotted 
black lines. We have explicitly checked the distribution of 
X 2 values for the individual mock catalogs to ascertain that 
the Gaussian approximation for the likelihood is valid, and 
found that the cumulative distribution of y 2 matches the ex- 
pected distribution at extremely high confidence (using the 
Kolmogorov-Smirnov test) at all points on the grid. Conse- 
quently, the method we have used to determine confidence 
regions is valid. The minimum A m i„, appears at f cen — 0.74 
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Figure 6. A, defined by Eq, lll6[ l and derived by comparing the projected 2PCCF ui p (r p ) for the whole sample of radio AGN as predicted 
by the halo occupation model and as measured from the SDSS data, is plotted with respect to the minimum value A m i„, in the grid of 
the two model parameters, f cm and M™ m . The contour levels, as indicated at the right-hand side, are chosen to cover the full range of 
A produced by all the models. The 68.3%, 90% and 95.4% confidence levels are plotted as solid, dashed and dotted black lines. 



and M^ in = 10 12 ' 75 h" 1 M Q with x'/d.o.f. = 0.5, indicating 
that the fit is acceptable. 

However, there is a strong degeneracy between the two 
parameters in the sense that models with smaller f cen and 
lower A/™ m can also provide a reasonable fit to the data. 
While this minimum is the preferred solution at the l-cr 
level, there is a banana-shaped degeneracy region extend- 
ing down to (0.60, 10.0) that is allowed at the 2-er level. The 
slight saddlepoint at (0.60, 11.75) that appears to divide this 
region does not have a sufficiently large AA relative to the 
minimum that we can robustly consider it is being real. As 
shown in Fig. [7] the residual noise in A due to the use of 
finite N moc k for the modeling is comparable to the size of 
AA that creates this apparent saddlepoint. The noise in A 
was determined by bootstrapping the (roughly independent) 
N m ock to make many new sets of mock catalogs and de- 



termining the variance between A for these sets of N moc k- 
However, the size of A A that distinguishes this degeneracy 
region from the rest of the plane in (fcen, M™ ln ) is signif- 
icantly larger than the noise in A. The reason that we are 
able to fit the data reasonably well at other points in this 
degeneracy region, despite the very different halo model pa- 
rameters, is that the higher satellite fraction places more 
AGN as satellites in very massive halos, offsetting the lower 
bias for central AGN due to the lower Af ™ in . 

We now use the center and right panels of Fig.f4]to eval- 
uate whether the model that best fits the data for the full 
sample can also describe the data for the sample split into 
two stellar mass bins. We have repeated the above analysis 
with two samples of radio AGN with logM, < 11.44 and 
logAf, ^ 11.44. We evaluate the correlation function using 
the same radial bins for 200 mock catalogues. We see that 
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Figure 7. Noise in contour plot of A, Figure [6] due to finite 
Nm OC k i a s determined from bootstrap resampling of the 200 mock 
catalogs. The noise as a function of N mock is shown for three 
different points on the grid: the global minimum, and two other 
points in the degeneracy region. 



while the model is able to describe the data for the lower 
stellar mass bin quite well, including the significant sepa- 
ration in signal between radio AGN and controls, there is 
some tension between the observations and the model for 
the higher stellar mass bin. As shown, the observed signals 
still have a significant offset, but the model signals are nearly 
the same for the radio AGN and for the controls. This is not 
surprising, since at these high stellar masses, essentially all 
halos are above M™ ln . Consequently, the apparent failure 
of the model at high stellar mass most likely results from 
the fact that a step- function model for the radio AGN prob- 
ability is overly simplistic. A probability that is a function 
of mass would allow for a better description, but unfortu- 
nately the data quality do not justify adding additional halo 
model parameters at this time, so we defer such an analysis 
to future work with more data. 



4-2.2 Galaxy-galaxy lensing 

Here we present the galaxy-galaxy weak lensing signal for 
the radio-loud AGN sample, along with the halo model fits. 
These results are shown in Fig. [S]for the full radio-loud AGN 
and control samples (upper left); a stellar mass subsample 
containing the lower 2/3 of the sample in stellar mass (lower 
left); the remaining upper 1/3 of the sample in stellar mass 
(upper right); and the radio- loud AGN split by log (P/Mbh) 
(lower right). 

A number of trends are evident in Fig. [5] First, we see 
that the lensing signal is higher on all scales (R < 2/i _1 Mpc) 
for the radio-loud AGN than for the control sample. Sec- 
ond, this trend persists in the different stellar mass subsam- 
ples. Third, the lensing signal increases with stellar mass, 
as expected. Finally, the lensing signal for the two bins in 



log (P/Mbh) are not markedly different from each other. For 
the lower and higher bins in this quantity, the mean stellar 
masses are 2.45 and 2.75 x 10 11 M@, respectively, a difference 
of only 10 per cent, much less than the difference between 
the mean stellar masses for the two stellar mass bins, so this 
near equivalence is not surprising. 

Next, we present the halo model interpretation of these 
results. Table[2]gives the best-fitting halo model parameters 
for each of the subsamples shown in Fig. [8] We focus on the 
results for central halo masses, because the satellite fractions 
are quite noisy. For the full sample, the best-fitting central 
halo mass M cent is 80 per cent higher for the radio-loud AGN 
than for the control sample at fixed stellar mass and redshift. 
When accounting for statistical correlations between the two 
samples using the bootstrap method, we find that the central 
radio-loud AGN have higher mass than the controls at the 97 
per cent confidence level. For the lower stellar mass bin, this 
conclusion is less statistically significant (50 per cent higher 



mass, with p(M t 



(AGN) 



> M, 



(controls) \ 



0.71) but for the 



higher stellar mass bin, it is more significant than for the 
full sample (factor of two higher mass, with p(M t 



(AGN) 



M^°" tro!s ') = 0.99). This result is in sharp contrast to the 
results for central optical AGN, which appear to have the 
same halo mass as the optical control galaxies. 

When comparing the amplitude of the central halo 
masses to tho se for the general galaxy population studied 
via lensing in iMandelbaum et al.l (|2006t ) , we find that the 
results for the control sample are consistent with that pa- 
per once we account for the different halo mass definitions, 
provided that we compare to the early type galaxy sample 
(which has higher mean central halo mass than the late type 
sample for M« > 10 n Af©). 

Given the high mean halo mass inferred for central 
radio-loud AGN, we must check our modeling assumption 
that the range of halo masses is narrow. Fig. shows a 
plot of the halo mass distribution for model B described 
in section 14.2.11 which is able to reproduce the radio AGN 
clustering results. As discussed, this HOD includes a mini- 
mum halo mass cutoff of 10 12,5 /i _1 Mq . As can be seen, the 
central halo mass distributions are indeed broader than the 
factor of ~ 6 that we need for the M cen t derived from lens- 
ing to be physically meaningful. As before, we do not apply 
a correction factor, due to the uncertainty in determining 
it. We merely note that since the AGN are more likely than 
the controls to reside in clusters (a claim that we will back 
up via direct comparison with a cluster catalogue), the cor- 
rection factor to get the mean central halo mass should be 
higher for the AGN than for the controls. Thus, the size 
of the central halo mass difference between radio-loud AGN 
and control galaxies is in fact underestimated by our neglect 
of these corrections. 

To avoid uncertainties in the best-fitting halo masses 
caused by the cluster membership of some galaxies, we have 
cross-correlated our radio-loud AGN and control samples 
with a pre-existing SDSS sample of galaxy clusters. We can 
then compare halo masses for radio-loud AGN and control 
galaxies in the field (i.e. excluding cluster members). For 
this sample, the distribution of halo masses will no longer 
have a tail extending to very high masses. 

Since the radio-loud AGN sample is predominantly (85 
per cent) in the redshift range 0.1 < z < 0.3, the natural 
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Figure 8. The galaxy-galaxy lensing signal for the radio-loud AGN and control galaxies split into subsamples as indicated on the plot. 
Points show the measured signal and lines show the best-fitting halo model. 



Table 2. Best-fitting halo model parameters for fits to the radio-loud AGN g-g weak lensing signal, with 68 per cent CL errors (in each 
case, marginalized over the other HOD parameter). 



Sample 



Radio- Loud AGN 



Me. 



[lO^/i-iJWe] 



Controls 



M c , 



[lO^h^M®] 



Full 

log(Af„/Af e ) < 11.44 
log(M„/M ) J; 11.44 
log(P/M bh ) < 15.7 
log (P/M bh ) Z 15.7 



1.6 ±0.4 
O.sj 
4.9" 1 



5 +0.4 
-0.5 
a +0.7 

-0.9 

1.8 ±0.4 

,1+0.4 
-0.9 



1.4P 



0.22 ±0.11 

,+0.30 
-0.16 
,+0.15 



0.3V 



0-011o.oi 
1 3+° 15 

40+ 40 
u ' w -0.15 



+0.12 
0.10 
+0.11 
0.09 
2.5 ±0.4 



0.91 
0.56 



0.13 ± 0.05 

-,+0.05 
-0.07 

0.16 ± 0.10 



0.10"' 



0.97 
0.71 
0.99 
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Figure 9. The distributions of halo mass in the mock catalogues that are able to reproduce the clustering signal for radio-loud AGN in 
the full sample and our two stellar mass bins. Results are shown for the full sample (black lines), central galaxies only (red lines), and 
satellites only (blue lines). 



choice of cluster cata logue is the SDSS MaxBCG catalogue 
l|Koester et al.ll2007bl lah. which contains clusters in the red- 
shift range 0.1 < z < 0.3 that are selected based on the 
existence of a red sequence. 

We select radio-loud AGN and control galaxies in this 
redshift range and check whether they are within l/t -1 Mpc 
(physical projected separation) of a cluster, and within 
Az = ±0.04 of the cluster BCG. We note that the choice 
of redshift separation is a factor of ten larger than the ve- 
locity dispersion of even the very largest galaxy clusters. 
The motivation for this choice comes from the typical pho- 
tometric redshift error of the galaxies in the maxBCG cat- 
alogue, and ensures that 95 per cent of true cluster mem- 
bers would be found. For reference, the minimum mass of 
the public maxBCG catalogue, with scaled ric hness ^ 10, 
is ~ 6 x 10 13 /i _1 M slln (defined using M 20 op, iReves et al.l 
2008). We also determined whether the radio- loud AGN or 
control is the cluster BCG or a satellite. These statistics are 
presented for radio-loud AGN and controls for both the full 
sample and for subsamples selected according to stellar mass 
in Table |3 

The numbers presented there for each sample are de- 
fined as follows: 

Number in sample that are in a cluster . 

/dust = r: r : ; (!') 

Number m sample 

and 

, Number in sample that are BCG of a cluster . . 

/bcg = r : i I 18 ) 

JN umber m sample 

Based on this table, we note a few interesting trends. 
First, for all subsamples we considered, /ciust and /bcg are 
higher for the radio-loud AGN than for the controls. This 
difference is most pronounced for the lower stellar mass bin, 
with a factor of two difference between radio-loud AGN and 
controls. For the higher stellar mass bin, the differences are 
at the ~ 20 per cent l evel. These resul ts are broadly in 
agreement with those of lBest et al. I (|2007h . 



Second, if we compare these two numbers for a given 
sample, we can determine the fraction of those in clusters 
that are BCGs, where the rest are satellites. For the full 
radio-loud AGN and the control samples, this number is 60 
per cent. Thus, on average, the distribution of centrals ver- 
sus satellites for those that are cluster members is the same 
for radio-loud AGN and control galaxies. For the lower and 
higher stellar mass samples, we again find consistency be- 
tween the radio-loud AGN and controls, with BCG fractions 
of those that are in clusters of 40 per cent and 85 per cent, 
respectively. 

It is apparent that the halo mass distributions of both 
the radio-loud AGN and control samples may have signifi- 
cant contributions from cluster BCG and satellite galaxies, 
which will skew the halo mass distribution to the high mass 
end. Consequently, we repeat the weak lensing analysis using 
only those galaxies in the redshift range 0.1 < z < 0.3 that 
are not within any maxBCG cluster. This cuts down the 
size of the sample significantly at the high stellar mass end, 
so we only analyze the full sample, without any divisions in 
stellar mass or radio power. 

The lensing signal for these "field samples" is shown in 
Fig. 1101 along with the best-fitting halo model. We initially 
used the same halo model as before, but found that the 
satellite fractions were all consistent with zero within the 
noise (as one would expect given the sample design). Thus, 
we redid the fits with fixed a = 0, fitting for central halo 
mass only. 

The mean stellar mass for these field samples is 2.1 x 
1O 11 M0. The reduction in best-fitting central halo mass go- 
ing from the full radio-loud AGN and control samples, to the 
field subsamples, is roughly 15 per cent. This reduction is 
expected, since we have excluded the AGN in clusters, which 
have typical halo masses of > 1O 14 M0. In both cases, how- 
ever, for the full samples and for the field subsamples, the 
best-fitting central halo mass for the radio-loud AGN sam- 
ple is roughly twice that of the control samples. We have 
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Table 3. Results of matching the radio-loud AGN and control samples with 0.1 < z < 0.3 against the maxBCG cluster catalogue. 
Quantities presented in the table are defined in Equations 1171 and 1181 



Sample 


/ oto8t (AGN) 


/dust (controls) 


/bcg (AGN) 


/bcg (controls) 


Full 


0.24 


0.16 


0.14 


0.09 


logM* < 11.44 


0.20 


0.12 


0.08 


0.05 


logitf* ^ 11.44 


0.36 


0.29 


0.31 


0.23 



Field galaxies only 




R [h-'Mpc] 



Figure 10. The galaxy-galaxy lensing signal for the radio-loud 
AGN and control galaxies located in the field with 0.1 < z < 0.3. 
Points show the measured signal and lines show the best-fitting 
halo model. 

already explored this result for the full sample; for the field 
sample, we find best-fitting central halo masses of (1.5±0.3) 
and (0.76 ± 0.14) x 10 l:i h' 1 M Q (radio AGN and controls, 
respectively). The difference between the two masses is thus 
significant at the 97 per cent CL. 

4.3 Joint constraints on halo masses and satellite 
fractions of radio AGN 

In section l4.2,ll we presented a best-fitting HOD model that 
well reproduced the clustering signal of radio AGN and the 
control sample. In this section, we compare the results of 
the model to the results obtained from the lensing analysis. 

In Fig. 111! we plot the mean halo mass of the central 
radio-loud AGN and the fraction of satellite AGN, as pre- 
dicted by the best-fitting model describing the clustering 
results, as a function of stellar mass. The results from the 
lensing analyses are also plotted. As discussed in the pre- 
vious section, the lensing analysis finds a higher mean halo 
mass for radio AGN at slightly more than 2a, which persists 
even when radio AGN residing in clusters are excluded from 
the analysis. This difference is even more significant, at the 
99 per cent CL, for the higher mass subsample. 

As shown, the model that best describes the clustering 
results is consistent with the g-g lensing results in the sense 



that different mean masses are predicted for the controls and 
the radio AGN. As already discussed, the fact that the same 
masses are predicted for the two samples in the higher mass 
bin is a consequence of overly-simplistic modeling. While the 
g-g lensing estimates of the satellite fraction are fairly noisy, 
they are also consistent with the model that describes the 
clustering. It would be valuable to confirm our conclusions 
using larger radio AGN samples that will be available in the 
future. This is particularly true at very low stellar masses 
(< 10 11 AIq), where a measurement with lensing was not at 
all possible and where we would be most sensitive to the dif- 
ference between radio AGN and controls. When samples are 
available with significantly better statistics, the data qual- 
ity will warrant a more careful analysis using the mocks to 
compare against both the clustering and lensing signal, and 
possibly more sophisticated halo modeling than that which 
was attempted here. 



5 SUMMARY 

We now present a comparison of the inferred halo masses 
and satellite fractions of both optical and radio- loud AGN, 
from both the clustering and g-g lensing techniques. 

In the left panel of Fig. 1121 we plot the inferred central 
halo masses of optical AGN (cyan) and radio AGN (green) 
from the lensing analysis as a function of stellar mass. The 
solid symbols give results for the AGN, while the open sym- 
bols give results for the control samples. 

In the middle panel, we plot the inferred central halo 
masses of optical (blue) and radio AGN (red) from the clus- 
tering analysis. The results are from the published model of 
iLi et alj l|2006bl ) for the optical AGN, and the magenta lines 
show our best-fitting model for the radio AGN presented in 
section 14.2.11 The halo mass as a function of stellar mass 
is plotted as solid lines for the AGN and as dashed lines 
for the control galaxies. As can be seen, the clustering tech- 
nique allows us to extend the halo mass measurements over 
a significantly larger range in stellar masses. We have ex- 
perimented with subsamples with even lower stellar masses 
and we found that we are able to derive a (noisy) clustering 
measurement for radio galaxies with stellar masses less than 
1O 11 M0, but there are too few objects to permit a lensing 
analysis to be carried out. Lensing measurements at lower 
stellar mass with larger, future datasets may be critical for 
reducing the modeling uncertainty indicated in this figure. 
Note that the convergence of the curves at high stellar mass 
is a consequence of our very simple model which has little 
or no flexibility to adjust the clustering of high mass stellar 
mass objects, the great majority of which are central galax- 
ies and live in massive halos. 

Finally, in the right-hand panel, we superimpose the 
clustering and lensing results. We see that the agreement 
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Figure 12. Left: Average central halo mass derived from galaxy-galaxy lensing is plotted as function of stellar mass, for three subsamples 
of optical AGN (cyan) and two subsamples of radio AGN (green). Filled symbols show results for the AGN samples, whole open symbols 
show results for the control samples. The stellar mass bin widths are indicated with horizontal errorbars. Middle: Average central halo 
mass derived from the clustering analysis is plotted as a function of stellar mass for optical (blue) and radio (magenta) AGN, and for the 
control galaxies (black dashed line for optical and red dashed for radio). Right: The results from the lensing and the clustering analyses 
are plotted on top of each other for comparison. 



between the results for these two completely independent 
techniques is satisfactory, once we take into account the fact 
that (a) there is some modeling uncertainty in the best- 
fitting central halo masses from lensing due to the assump- 
tion of a narrow central halo mass distribution, and (b) the 
halo masses from the lensing analysis are relatively indepen- 
dent of the assumed erg (which to lowest order only affects 
the best-fit a) , but the halo masses from the clustering anal- 
ysis are tied to as — 0.9 from the Millennium simulation 
through the large-scale bias - halo mass connection. In the 
latter case, we can estimate the effect of lowering as from 
0.9, as in the mocks, to 0.8 (as in the WMAP 5-year results, 

0.25. Using the mass 



iDunklev et all 



2008) at fixed Q r , 



Tink er et all (|2008f ) at the typical redshift 



function from 

of radio AGN, and translating to our halo mass definition, 
we find that the requirement that the stellar mass function 
be matched, which is essentially an abundance constraint, 
would lead to the masses from the clustering analyses be 
lowered by ~ 15 per cent at 10 13 h~ x Mq, or ~ 25 per cent 
at 3 x 1O 13 /i _1 M . While this difference is significant, it can 
only partially account for the differences shown in the fig- 
ures. However, the modeling uncertainty due to the simple 
HOD used for the lensing analysis can lead to significant 



additional uncertainty in those masses, typically leading to 
underestimation (i.e., the sign of the apparent discrepancy) 
by several tens of per cent when the central halo ma ss dis- 
tributions are quite broad (|Mandelbaum et~a l. 2005b), as is 
the case for several of our stellar mass subsamples (Figs. [3] 
and©. 

Putting the lensing and the clustering results together 
leads us to the following major conclusions: 

• Radio AGN are hosted by galaxies with higher stel- 
lar masses than optical AGN, and are also in more mas- 
sive dark matter halos. The mean stellar mass of the optical 
AGN sample is 8 x 1O 1O M0 and the corresponding mean 
central halo mass deduced from galaxy-galaxy lensing is 
(8.0T1.5) x 10 11 /i" 1 M Q . The mean stellar mass of the radio 
AGN sample is 2.5 x 10 1 Mq and the mean halo mass de- 
duced from galaxy-galaxy lensing is (1.6±0.4) x 10 13 hT 1 Mq. 
Thus, the mean stellar mass and mean central halo mass for 
the radio AGN are ~ 3 and ~ 20 times the corresponding 
values for the optical AGN. 

• At fixed stellar mass, radio-loud AGN inhabit more 
massive dark matter halos than optical AGN. This is seen 
both in the clustering and and in the galaxy-galaxy lensing 
analyses. Note that for the g-g lensing analysis, the high- 
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Figure 11. Mean halo mass for central radio AGN (upper panels) 
and satellite fraction of all radio AGN (lower panels) are plotted 
as function of stellar mass. In each panel, the red thick line shows 
the best-fitting model determined by the clustering measurements 
of radio AGN, and the blue thick line shows the result for control 
galaxies. The thin lines show the 1 — a variance between 200 mock 
catalogues. The results determined by the g-g lensing analyses are 
plotted as magenta circles for radio AGN and as green triangles 
for control galaxies; horizontal errorbars indicate the widths of 
the stellar mass bins. 



est stellar mass bin for the optical AGN and the lowest bin 
for the radio-loud AGN have the same mean stellar mass. 
Fig. [12] shows that these two samples have very different 
halo masses. While the halo masses are also different for the 
control galaxies due to the morphology- dependence of halo 
mass at fixed stellar mass > 1O 11 M0 (|Mandelbaum et al.l 



2006), the difference is even more pronounced for the opti- 
cal and radio-loud AGN samples. 

• At fixed stellar mass, optical AGN inhabit dark matter 
halos of similar mass as galaxies of the same stellar mass 
selected without regard to AGN properties. 

• At fixed stellar mass, radio-loud AGN inhabit more 
massive dark matter halos than galaxies of the same stellar 
mass selected without regard to AGN properties. We em- 
phasize that despite the difficulty in representing this offset 
using our simple clustering model, it is present observation- 
ally at high significance in both the clustering data (Fig. [4| 
and the lensing data (Figs. [5] and 1 10 p both for the full sam- 
ples and for the two stellar mass subsamples. 

• The clustering and lensing analyses together favour a 
model in which radio-loud AGN are not found in dark matter 
halos with masses less than about 3 x W 12 h -1 Mq, though 
the preference for this over a model with a minimum mass 
a factor of 10 smaller is only a 2a difference given the size 
of current datasets. 

One unresolved point in the reconciliation between the 
modeling of the clustering and its comparison with the lens- 
ing analysis is that the lensing analysis found the most sig- 
nificant difference between the central halo masses for radio 
AGN and controls at the high stellar mass end, whereas the 
clustering modeling suggests the largest difference should 
occur for lower stellar mass. This point is simply an artifact 
of overly simplistic modeling of a strict mass threshold, and 
may therefore be resolved with more sophisticated modeling 
involving a probability that is a function of mass; however, 
the introduction of more model parameters is not justified 
by the achievable S/N of the data at this time. 

In Fig. 1131 we compare the satellite fractions inferred by 
the two techniques. The results are broadly consistent with 
each other. The clustering analysis yields bigger differences 
in the inferred satellite fractions, particularly for AGN with 
low stellar masses. Optical AGN with low stellar masses are 
predicted to be quite strongly biased to central galaxy hosts, 
whereas radio AGN with low stellar masses are predicted to 
be located more frequently in satellite galaxies. The lensing 
analysis gives some weak indications of trends in the same 
direction, but the results are far from conclusive. 



6 IMPLICATIONS OF THIS WORK 

Perhaps the most important finding of this work is that op- 
tical AGN largely follow the same relation between stellar 
mass and halo mass as "ordinary" galaxies, but that radio- 
loud AGN deviate significantly from it (by a factor of ~ 2). 
This statement is true at the 97 per cent CL even if we 
restrict the analysis to those radio-loud AGN and control 
galaxies in the field (excluding those in massive groups and 
clusters) . 

This result implies that the large-scale halo environ- 
ment plays an important role in under standing the radio 
AGN phenomenon. Previous work (e.g. iBest et all [2005a) 
has shown that the fraction of radio AGN increases strongly 
for more massive galaxies, suggesting that radio jets are 
more readily triggered in galaxies with more massive black 
holes. However, when we compare radio AGN with control 
galaxies of the same stellar mass selected without regard 
for nuclear activity, we find that the radio AGN reside in 
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Figure 13. Satellite fraction as a function of stellar mass. Sym- 
bols and lines are the same as in the previous figure. 



dark matter halos that are a factor of two more massive 
on average. This boost in halo mass appears to be largely 
independent of luminosity of the radio source. This demon- 
strates that black hole mass is not the only parameter that 
controls the radio AGN phenomenon - some aspect of the 
larger-scale environment of the galaxy must play a crucial 
role in regulating when the jet is switched on or when it is 
visible at radio wavelengths. 

Recent semi-analytic models have assumed that feed- 
back from radio AGN only becomes important in halos in 
which gas is cooling quasi- statically, i.e. halos above a mass 
of a f ew xlO 11 - IO^Mq (|Croton et all 120061 ; iBower et all 
120061 : rCattaneo et al.ll2006h . Our clustering results strongly 
support this idea. Our two best-fitting HOD models both 
invoke a minimum halo mass close to these values, below 
which radio AGN are no longer found. If we do not impose 
a minimum mass, our models are not able to fit the ob- 
servation correlation amplitude of radio-loud AGN on large 
scales. 

Finally, the fact that the optical AGN in our sample 
follow the same M» — Mhaio relation as the general galaxy 
population, implies that the optical AGN phenomenon is 
largely decoupled from the host halo. For the highest stel- 
lar mass bin, the best-fitting central halo mass derived from 
galaxy-galaxy lensing is more consistent with published re- 
sults for late-type galaxies. This result is consistent with the 
general tendency of narrow-line AGN to be associated with 
galaxies with ongoing star formation and by extension, a 
cold interstellar medium. 
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